Time Reversibility of Quantum Diffusion in Small-world Networks 
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We study the time-reversal dynamics of a tight-binding electron in the Watts-Strogatz (WS) 
small-world networks. The localized initial wave packet at time t = diffuses as time proceeds until 
the time-reversal operation, together with the momentum perturbation of the strength n, is made 
at the reversal time T. The time irreversibility is measured by I = \H(t = 2T) — II(i = 0)|, where 
fl is the participation ratio gauging the extendedness of the wavefunction and for convenience, t is 
measured forward even after the time reversal . When n = 0, the time evolution after T makes the 
wavefunction at t = 2T identical to the one at t — 0, and we find 1 = 0, implying a null irreversibility 
or a complete reversibility. On the other hand, as rj is increased from zero, the reversibility becomes 
weaker, and we observe enhancement of the irreversibility. We find that I linearly increases with 
increasing n in the weakly-perturbed region, and that the irreversibility is much stronger in the WS 
network than in the local regular network. 
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Time- irreversible phenomena in nature, such as spread 
of an ink blob in water and aging of living organisms, are 
ubiquitous in our macroscopic world We always see 
the ink blob spread in a teacup and people get older, not 
the other way around. These time-irreversible behaviors 
look puzzling because all these dynamics are based on mi- 
croscopic equations of motion, quantum or classical, that 
contain time-reversal symmetry. In the late 19th century, 
this seemingly paradoxical observation of time-reversal 
microscopic dynamics and time-irreversible macroscopic 
thermodynamics perplexed many scientists until a break- 
through was made by Boltzmann. It is now well-known 
that the directionality of the arrow of time in the second 
law of thermodynamics can be understood by the huge 
difference in the numbers of allowed microscopic states 
between the initial and the final macroscopic states. 

Recently, the quantum and the classical time-reversal 
dynamics of a diffusion system and a spin system have 
been studied to answer the long-standing question of 
where does the irreversibility come from In Ref . |, 

quantum diffusion systems were shown to display the uni- 
versal behavior of time irreversibility in terms of per- 
turbation strength. Because the complex networks have 
drawn much interest, there have been studies of quan- 
tum systems in Watts-Strogatz (WS) small-world net- 
works [6|. Localization phenomena of quantum systems 
in small- world networks have been studied 0, i|, and 
quantum diffusion problems have been investigated in 
various ways in WS networks [§-[l2j]. 

In the present work, we aim to study the time re- 
versibility of the quantum diffusion problem in WS net- 
works. The diffusion of a wave packet has been shown 
to occur much faster in WS networks than in a regular 
network [l(| [H| : A suitably defined diffusion time in a 



quantum system shows r ~ log N in the former while 
t ~ N for the latter structure. An abrupt change in the 
scaling behavior of r has been shown to occur as soon 
as the rewiring probability becomes nonzero. In other 
words, a transition from the slow world (r ~ N) to the 
fast world (r ~ log N) occurs at null rewiring probabil- 
ity simultaneously with a structural transition from the 
large world (I ~ N) to the small world (I ~ log N), with I 
being the average path length connecting two arbitrarily 
chosen vertices. Consequently, a change in the scaling be- 
havior of t reflects a dynamical aspect of the small- world 
transition in WS networks. In this work, we examine the 
effect of shortcuts on the time reversal dynamics of the 
tight-binding electron system at different perturbation 
strengths. 

The time evolution of the tight-binding electron in 
a network structure is governed by the following time- 
dependent Schrodinger equation: 
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where h is the Planck constant, and the quantum me- 
chanical kct \^>) and the Hamiltonian H in position rep- 
resentation are written as ^„ = (n\^) and 
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otherwise, 
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respectively. Here, n is the vertex index, A„ is the set of 
directly connected vertices of n, and the on-site energy 
has been assumed to be uniform and set equal to zero. As 
a further simplification, we also assume that the hopping 
energy A does not depend on n or n'. Henceforth, we use 
the dimensionless units: the time is measured in units of 
fi/A, the position in units of the lattice spacing a of the 
one-dimensional (ID) regular lattice without shortcut, 
and the momentum in units of h/a. 

Once the WS network for a given rewiring probability 
p is constructed following the standard procedure (see 
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Ref.Q), we perform the numerical integration of the time- 
dependent Schrodinger equation starting from the initial 
localized wave packet at time t = (*„ = S n ,o) by us- 
ing the fourth-order Runge-Kutta algorithm with the dis- 
crete time step dt = 0.01. Except for p — 0, all results 
are obtained from the average of 1000 different network 
realizations. We use the time-reversal test similarly to 
Ref. S At the reversal time T, a momentum perturba- 
tion of the strength rj is made by applying the operator 
exp(i?7x) to \^{T)) with the position operator x. For 
convenience, we measure t in the forward direction even 
after the reversal time T so that at t = 2T, the system 
goes back to the initial state in the absence of the per- 
turbation (i.e., when r) = 0). During the numerical time 
integration, we compute the participation ratio 
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which has been frequently used in the study of localiza- 
tion phenomena. When the wavefunction is completely 
localized, II takes the value of unity. On the other hand, 
as the wavefunction spreads in space and the quantum 
system becomes extended, II is known to be O(N). The 
time irreversibility I is then defined as 



i= |n(2T) -n(o)|. 



(4) 



It is to be noted that H(t < T) is independent of 77 
whereas n(2T) and I are functions of rj. When rj = 0, I 
takes the value zero, meaning the null irreversibility (or 
the complete reversibility), because U(2T) — 11(0). As 
77 is increased from zero, the perturbed momentum at T 
makes it difficult for the system to go back to the initial 
quantum state; thus, I is expected to increase. 

Figure Q] displays II(i) in Eq. © (a) for the local ID 
regular network corresponding to the WS network at the 
rewiring probability p — and (b) for the WS network 
at p = 0.2. Although not clearly discernible, each of 
Fig. [TJa) and (b) has 101 different curves corresponding 
to r) = 0.0,0.001,0.002, •• • ,0.1 (from bottom to top). 
We also show the time irreversibility I at various values 
of p in Fig. HJc). The system size N = 1600 and the 
reversal time T = 2 are used. We first observe that H(T) 
is about four times larger for p — 0.2 than for p = 0.0, 
which is explained by the fast diffusion in the WS net- 
work structure, as discussed in Ref. [To|- At both (a) p = 
and (b) p = 0.2, time reversal dynamics is shown to give 
rise to a larger reflection asymmetry of H(t) around t = T 
as i] is increased. This is not surprising because a larger 
momentum perturbation will surely make the time re- 
versed quantum state at t = 2T more different from the 
initial state at t = 0. A more interesting observation 
one can make from the comparison of Fig. [1] (a) and (b) 
is that the deviation of T1(2T) from 11(0) is much big- 
ger for the WS network (p — 0.2) than for the regular 
network (p — 0.0). This clearly indicates that the struc- 
tural irregularity in the WS network plays an important 
role in making the perturbed dynamics deviate severely 
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FIG. 1. (Color online) The participation ratio II (t) as a 
function of time t in (a) the local regular network and (b) 
the WS network at rewiring probability p = 0.2. The sys- 
tem size N — 1600 and the reversal time T = 2 are used. 
Although not clearly discernible, (a) and (b) contain 101 
different curves for the momentum perturbation strength 
■q = 0.0, 0.001, 0.002, ■ ■ ■ ,0.1 (from bottom to top). The per- 
turbed and time-reversed quantum state in the WS network 
in (b) shows a much stronger deviation from the initial local- 
ized quantum state, (c) The irreversibility i=|II(2T) — II(0)| 
versus rj at various rewiring probabilities p. When p = 0.0 
for the local regular network, / increases with increasing rj. 
On the contrary, when p 7^ 0.0, I first increases, but soon 
saturates. The oscillatory behavior observed for p — 0.0 orig- 
inates from the discreteness of the underlying periodic lattice 
structure. 



from the unperturbed one. The enhancement of the ir- 
reversibility in the WS network is more clearly seen in 
Fig. [TJ (c) for / versus 77 at various rewiring probabil- 
ities p. Figure [IJc) shows that the irreversibility is a 
monotonically increasing function of the rewiring prob- 
ability p. In other words, the enhanced structural ran- 
domness caused by more shortcuts drives the quantum 
diffusion dynamics in the WS network to become more 
sensitive to the momentum perturbation. Interestingly, 
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FIG. 2. (Color online) The time irreversibility / in Eq. ([4]) is 
measured for various values of the rewiring probability p at 
a fixed value of the reversal time T to plot I versus pT for 
(a) T = 1 and (b) T = 10. For pT < 1, the system is in 
the ID region and / does not depend on the system size N. 
As pT is increased to become pT > 1, the system enters the 
small- world region and different sizes begin to show deviations 
from each other. As pT is increased further, I saturates to 
some finite value: I(pT — > oo)/N ~ const., independent of 
T, as shown in (c). The strength of momentum perturbation 
is rj = 0.1, and all results are averages over 1000 different 
network realizations. 

the WS network at a sufficiently larger rewiring probabil- 
ity behaves very differently from the local regular lattice: 
For p > 0, I increases linearly with increasing 77 for the 
weakly-perturbed region and then soon saturates to a fi- 
nite value. In contrast, for p = 0, / increases linearly 
with increasing rj without showing saturation. However, 
these findings need to be carefully examined in view of 
the finite scales in the system, i.e., finiteness of the length 
scale (N) and the time scale (T). 

It should be noted that in WS networks there is an 
additional length scale other than the microscopic length 
scale of the lattice spacing a and the macroscopic length 
scale of the system size N. When a nonzero rewiring 
probability p is given, the number of shortcuts is pro- 
portional to pN , which determines the third length scale 



of the average distance £ between the two endpoints of 
shortcuts, i.e., £ ~ N/pN — l/p. If the tail part of the 
wavefunction of the tight-binding electron does not have 
enough time to arrive at the closest shortcut endpoint, 
the system should behave just like a ID regular lattice. 
Also, quantum diffusion in a regular ID system is known 
to have a diffusion time proportional to the system size to 
become fully extended (r ~ N) , which indicates that the 
wavepacket in the ID system spreads on a distance scale 
i ~ T before reversal time. Accordingly, if T ~ I < l/p, 
the system behaves just like a ID regular lattice while 
it changes its behavior as T ~ l/p is crossed. Let us 
consider several cases: (i) For pT < 1, the system be- 
haves like a ID system because T is short enough that 
the wavepacket does not have time to arrive at a shortcut 
endpoint. Furthermore, in this case of short-time diffu- 
sion, all observed results should be independent of N as 
long as T <C N; i.e., T is too short for the particle to 
feel the finiteness of the system, (ii) For pT > 1, the 
tight-binding particle now begins to meet shortcut end- 
points and the diffusion behavior changes. As soon as 
the WS network begins to have a finite fraction of short- 
cuts, the quantum diffusion behavior is known to rapidly 
change so that the diffusion time r increases only log- 
arithmically with increasing N. This implies that the 
wavefunction spreads in a distance that increases expo- 
nentially with time. We, thus, conclude that this very 
fast diffusion occurs when pT > 1 and that the parti- 
cle can arrive at the other side of the system, giving rise 
to an iV-dependence in the irreversibility. Consequently, 
the size-independence of / in the region of pT < 1 is 
changed as we enter the intermediate region of pT > 1 . 
(iii) For the limiting case of pT 3> 1, the wavefunction of 
the particle has already completely spread; thus, the ir- 
reversibility should not depend on pT. In this long-time 
limit, the wavefunction becomes fully extended, and the 
participation ratio is IT ~ N. Consequently, for a big 
enough system, I(pT)/N saturates as pT becomes larger, 
approaching a value that does not depend on N . 

The results for the irreversibility I versus pT curve at 
rj = 0.1 shown in Fig. [2] clearly fit the above hand- waving 
scaling arguments well. In Figs.[^a) and (b), / is shown 
to be independent of N for pT < 1 whereas we see a 
clear size dependence for pT > 1. As pT is increased 
even further, / saturates to some values that depend on 
N, as clearly seen in Fig. [5{b). Figure [^c) shows that 
the saturated value of I/N in the limit of pT 3> 1 is 
independent of N. 

In summary, we studied the time reversal dynamics of 
the tight-binding electron system in the WS small-world 
networks. Initially, the localized quantum mechanical 
state of the electron evolves in time until the momentum 
perturbation of the strength rj is made at T, at which 
time reversal operation is also made. The ID regular 
lattice and the WS small-world network exhibit different 
behaviors: the irreversibility, measured by the difference 
between the participation ratio of the final and the initial 
states, is found to be bigger in the latter network due to 
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the structural randomness caused by the random short- 
cuts. It is also found and argued that the irreversibility 
does not depend on the system size for pT < 1 and that 
it saturates to a value as pT becomes larger. 
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